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Abstract 



The proper support of generalized parton distributions from relativistic constituent 
quark models with pointlike constituents is studied. The correct support is guaran- 
teed when the vertex function does not depend on the relative minus-momentum. 
We show that including quark interactions in models with pointlike constituent 
quarks might lead to a support problem. A computation of the magnitude of the 
support problem in the Bonn relativistic constituent quark model is presented. 
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Historically, form factors and parton distribution functions have played a key 
role in hadron physics. In 1997, a new research tool emerged in the framework 
of the generalized parton distributions (GPDs) [T], providing us with a natural 
unification of form factors and parton distribution functions |2] . GPDs describe 
the soft, i.e. non-perturbative hadronic, scattering part of Hard Exclusive 
Meson Production (HEMP) and Deeply Virtual Compton Scattering (DVCS). 
Today, a lot of experimental and theoretical effort is devoted to the study 
of these generalized parton distributions, of which is hoped that they will 
shed light on the origin of the spin of hadrons in general, and the nucleon in 
particular. 

At present, calculating GPDs from first QCD principles is not feasible, and one 
relies on chiral perturbation theory [5], lattice QCD [1], parametrizations [S] 
or phenomenological models [51l7f51l51[TI]|llfl2lll3lll4|15|16j . This letter focuses 
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on the latter class. In order to compare the GPDs computed in a phenomeno- 
logical model with the DVCS and HEMP data, a Q^-evolution needs to be 
performed [TTj. The evolution equations depend highly on the kinematic re- 
gion of the process. More specifically, one can distinguish between the DGLAP 
and ERBL regions, where the Q^-evolution is respectively governed by the 
DGLAP and the ERBL equations. As indicated schematically in Fig. [H par- 
ton model constraints make GPDs vanish outside these two regions [T81IT9] . A 
model is said to have the correct support, when it can resolve the DGLAP and 
ERBL regions, and has a vanishing GPD outside this region. Developing such 
a model is far from straightforward. In this letter we will show how the sup- 
port can arise naturally in quark models that are based on the Bethe-Salpeter 
approach. Including quark dynamics in the model, however, is likely to destroy 
the correct support. 

In the following, we will focus on the GPD of a pseudoscalar meson such 
as the pion, considering the case where the partons do not transfer helicity. 
Pseudoscalar mesons have only one such generalized quark distribution asso- 
ciated with them, whereas the nucleon has four. The theoretical analysis of 
the support region is, therefore, less cumbersome for the pion GPD than for 
the nucleon ones. 

GPDs are non-diagonal matrix elements of a bilocal field operator on the light 
cone. For partons with spin 1/2 (quarks) in a pseudoscalar meson, the GPD 
associated with helicity conserving partons is defined as follows [IS]: 

In this equation, P = ^ . Definition ([T]) uses light-cone coordinates. A 
fourmomentum p = [p^ , , , p^) can be written in light-cone coordinates as 
p = [p'^ ,p~ ,p_i] with p"^ = (p° ±p^)/\^ and = 

In order to calculate the GPD, several choices regarding the kinematics of the 
problem have to be made. In fact, two choices are used in literature [20|2T] . 
We make use of the definitions from [20] , illustrated in Fig. [2] for the DVCS 
process: x denotes the fraction of the average pion plus-momentum that is 
reabsorbed by the meson, while the skewedness ^ is a measure for the plus- 
momentum that is lost in the process: 

p+ _ p>+ 

^ ~ p+ _|_ p'+ ' '-^'^ 
^ is defined on the interval [0, ^max] with 




where t = = (P' - P)^. The abovementioned DGLAP and ERBL regions 
are defined as follows: 

• X G [—1, — DGLAP region (emission and absorption of antiquark); 

• X G [— ^,^]: ERBL region (emission of both quark and antiquark); 

• X G 1]: DGLAP region (emission and absorption of quark). 

This is depicted schematically in Fig. [T] As mentioned before, the GPD should 
vanish outside these regions, i.e. for |x| > 1. 

In a Bethe-Salpeter quark model, the bilocal current matrix element from 
Eq. ([T]) can be written in terms of the Bethe-Salpeter amplitudes xp'- 

= I d% tTr {xpK4, |) ((^ ^4 - ^i)Xp{x2, -|)) 7+} (4) 
+ Jd% zTr{((z ^y'^-m,)xM-"-,y[))l^XpC^,y[)] • 

The first (second) term on the right-hand side of this equation refers to the 
coupling of the bilocal current to the antiquark (quark). Both terms can be 
treated in exactly the same manner. In the remainder of this letter, we focus on 
the analysis of the second term of Eq. (j4]). Inserting the quark term in Eq. ([1]) 
and making use of the definition of the quark and antiquark propagators, 

{z ^^-m)S^{x) = i6^'^\x-x'), (5) 

we can rewrite the quark GPD as an integral in momentum space: 

""-^""'^'^^ 2 7 (27r)4^ 1^ 2(1 + ; 

X Tr |f p,(p + j)S[{y +P+ + P)^p(P)S2i-^ + P)} • 

(6) 

On the left-hand side of Eq. ([HD, the superscript q indicates that only the 
coupling to the quark is considered. The vertex functions 

Tp = {S[)-'XP{S^)-' (7) 

have been introduced on the right-hand side of the equation, with subscript 1 
(2) referring to the quark (antiquark). The pictorial representation of Eq. ([6]) is 
sketched in Fig.[3l In some Bethe-Salpeter models, see e.g. Ref. [22], the quark 
and antiquark propagators are assumed to have the form of a free propagator 
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of pointlike fermions, with a constituent mass m: 



Sl,2iP) = 2 ^ ■ • (8) 

As it appears, the denominators of the three propagator terms in Eq. ([6]) 
suggest immediately the correct support region for the GPDs. Since the poles 
of the propagators are Lorentz invariants, any convenient frame can be chosen 
to prove this statement. The GPD in Eq. ([T]) was defined for a coordinate 

system with P along the z-axis [1]. The Breit frame is such a coordinate 
system, with P = (M, -|) and P' = (M, |). Here, = - | and A = 

(A_L, — 2^M) with ^ as introduced before. We emphasize that these expressions 
uniquely apply to the calculation of the elastic GPD, M' = M. 

Making use of the 5-function to perform the integration over p"*", it is now 
possible to write the propagator denominators as: 



P 
~2 



- + vf - + ze) = i3 + ^M[x + i)v~ + ie , 
- Y ^ - ml + iej = C + V2M{x - l)p~ + ie , 

where A, B and C do not depend on p~ . From these expressions, it is clear 
that the propagators have poles in the complex p~-plane. The position of these 
poles with respect to the real p~-axis (upper or lower half-plane) depends on 
the kinematics, more specifically on the sign of (x — ^), (x + ^) and (x — 1). 
For X ^ [—^,1], the propagator poles lie on the same side of the real axis. 
Furthermore, the propagators are analytic in all other points of the complex 
plane. Assuming that the vertex function P does not contain poles in the 
complex p~-plane, and making use of Cauchy's theorem, one can conclude that 
the momentum dependence of the propagators ensures the correct support for 
the GPDs. 



The above considerations are clearly valid if the vertex function P and its 
adjoint do not depend on p~ and are therefore free from poles when analy- 
tically continued in the complex p~-plane. The calculation of the pion GPD 
with constant vertex functions has been performed by S. Noguera et al. in the 
NJL-model [TU] with a Pauli-Villars regularization procedure. The calculation 
of GPDs with ^"-independent vertex functions was performed by H.-M. Choi 
et al. [T21IT3] . as well as by B.C. Tiburzi and G.A. Miller [T1IT51IT6] . These 
authors made use of a reduction of Bethe-Salpeter amphtudes to light-cone 
wave functions by projection on the light-cone. In Refs. [ri|15fl6j . the scalar 
Wick-Cutkosky model was adopted and successfully applied to the calculation 
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of scalar meson GPDs. In Refs. [T21IT3] . the light-front Bethe-Salpeter vertex 
functions were replaced by wave functions obtained in a light-front constituent 
quark model. All of these calculations yielded the correct support, in agree- 
ment with the present analysis. 

In instant-form models, however, the vertex functions can depend on p~ , which 
makes them have a different pole structure. In this case, the correct support 
can no longer be guaranteed. To prove this, we start with Liouville's theorem, 
which states that the only bounded entire functions are the constant functions. 
As a result, a vertex function which depends on the momentum variable p~ 
will not be bounded for all points in the complex p~-plane. Moreover, the 
radial part of a vertexfunction for a ground-state particle is a real-valued 
function on the real p~ axis, up to a constant phase factor. This is an important 
constraint. Indeed, for a holomorphic function / whose restriction to the real 
numbers is real-valued it can be proven that f{p~*) = f*{p~)- Here, * stands 
for complex conjugation. In other words: if / has a singularity in the upper 
half-plane, it also has one in the lower half-plane, and vice versa. Combining 
these two statements, we find that a dynamical, momentum dependent vertex 
function will either have complex conjugated poles, or a singularity on both 
complex half-circles at a; ± zoo. In the latter case, Cauchy's theorem is no 
longer applicable, while in the former case, there will always be at least one 
pole in each half-plane. Consequently, it is not a priori clear whether the 
GPD will vanish outside the interval x G [—1, 1]. This means that the correct 
support can no longer be guaranteed. Similar arguments apply to relativistic 
quark-diquark models of the nucleon with pointlike constituents. 

We stress that this does not mean that there will necessarily be a support prob- 
lem in Bethe-Salpeter constituent quark models with pointlike constituents 
containing dynamics - it merely shows that a support problem can no longer 
be excluded. The above considerations indicate that guaranteeing the correct 
support puts non-trivial constraints on the analytic properties of the Bethe- 
Salpeter vertex functions. Whether these are fulfilled in particular models has 
to be checked in each case individually. As an example, we show in Fig. Hlthe 
pion GPD calculated in the framework of the covariant quark model of the 
Bonn group [22] . 

In the Bonn model, the interaction kernels are instantaneous, i.e. independent 
from relative energy variables in the meson rest frame. More specifically, the 
employed interquark interactions are the linearly rising confinement interac- 
tion and the 't Hooft instanton induced interaction [23]. The confinement in- 
teraction has an appropriate Dirac structure to optimize the splitting between 
spin-orbit partners. The 't Hooft instanton induced interaction acts as a resid- 
ual interaction and accounts for the mass splittings in the pseudoscalar and 
scalar sectors [21]. Together with the assumption of free fermion propagators 
with an effective mass (Eq. (jHD), the instantaneous approximation allows to 
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reduce the full (but unsolvable) Bethe-Salpeter equation to a solvable Salpeter 
equation [25]. The resulting vertex function is independent of the relative en- 
ergy variable, but does depend on the relative minus- momentum. The Bonn 
model has only seven parameters which are fitted to the static properties of 
the meson spectrum. In Refs. [22|25|26|l27] . the model is described in detail. 

Figure H] shows the result of a calculation of the pion up-quark GPD if "+ as a 
function of x at t = —0.1 GeV^ and three different values for ^ {C, = 0, C, = 0.3 
and = 0.7). The first moment of the GPD, 



oo 



dxH^+{x,^,t) (9) 
yields the electromagnetic form factor, while the isospin symmetry relation 



H:+{x,U) = -HtA-x,^,t), (10) 

is exactly fulfilled. Here, H!^+ is calculated from the first term on the right- 
hand side of Eq. (jlj), which describes the coupling of the bilocal current to 
the antiquark. As a measure of the support problem, the support parameter 
(j) can be introduced: 

When the appropriate kinematical regions are resolved (see Fig. [1]), this frac- 
tion equals one. In contrast, (p < 1 implies that the GPD is non-zero outside 
the supported region, which is [— ^, 1] for the quark term. The smaller the frac- 
tion, the worse the support. In the case of Fig. HJ this fraction equals (j) = 0.13 
for ^ = and C, = 0.3. For ^ = 0.7, one has (p = 0.12. From these numbers, 
it is clear that the correct support is violated in the Bonn constituent quark 
model. 
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DGLAP ERBL DGLAP 

^ ^ ^ ^ ^ X 

-1-^0^ 1 

Fig. 1. Pictorial indication of the supported regions. 




Fig. 2. Feynman diagram of DVCS, with definition of the kinematical variables x, 
^ and t. 




Fig. 3. Schematic representation of the expression of the pion quark GPD in terms 
of Bethe-Salpeter vertex functions, Eq. ([6]). 



0.6r 




Fig. 4. The pion GPD H^^ as a function of x, calculated in the Bonn constituent 
quark model at t = —0.1 GeV^. Results are shown for ^ are ^ = (solid line), 
^ = 0.3 (dashed line) and ^ = 0.7 (dotted line). 
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